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Abstract 

Molecular dynamics simulations are used to study the mechanism and 
kinetics of hydrated electron diffusion. The electron centre of mass is 
found to exhibit Brownian-like behaviour with a diffusion coefficient con- 
siderably greater than that of the solvent. As previously postulated by 
both experimental and theoretical works the instantaneous response of the 
electron to the librational motions of surrounding water molecules consti- 
tutes the principal mode of motion. The diffusive mechanism can be un- 
derstood within the traditional framework of transfer diffusion processes, 
where the diffusive step is akin to the exchange of an extra-molecular 
electron between neighbouring water molecules. This is a second-order 
process with a computed rate constant of 5.0 ps _1 at 298 K. In agreement 
with experiment the electron diffusion exhibits Arrhenius behaviour over 
the temperature range 298 K-400 K. We compute an activation energy of 
8.9 kj moP 1 . Through analysis of Arrhenius plo ts and the application of 
a simple random walk model it is demonstrated that the computed rate 
constant for exchange of an excess electron is indeed the phenomenological 
rate constant associated with the diffusive process. 
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1 Introduction 

The hydrated electron has long been a topic of fundamental interest [H HI El E] ■ 
Since its original identification [5] it has become a focus of study particularly 
pertinent to radiolytic and electron initiated processes [5] . Ionisation of a polar 
solvent, by radiolysis or photolysis, leads to the production of a solvated electron. 
This transient species may localize and thermalize in the solvent, eventually 
being captured in a potential energy well formed by the solvating molecules. 
Whilst being all the time solvated, fluctuations in the solvent coordinate enable 
diffusion through the solvent medium. Upon encountering an acceptor species 
of sufficient electron affinity the electron may be incorporated into one of its 
orbitals. In such cases the rate of electron transfer or electron-initiated processes 
is inextricably linked to the diffusive dynamics of the solvated electron [7]. 

The hydrated electron, identified by its characteristic optical absorption 
spectrum [5, 8 , has been closely studied in ambient conditions by pump-probe 
spectroscopy [51 [TUl HH H3 HH ESI HI! HZ1- These works reveal a sophis- 
ticated picture of the early-time, ultra-fast dynamics of solvation from local- 
isation to equilibration and have been well complemented by computer simu- 
lations [13 H3 [2D]. The diffusive dynamics following equilibration have also 
been a focus of interest. Conductometric measurements of electron diffusion in 
water reveal both an unusually high mobility (several times greater than that 
of classical ions) and Arrhenius-like diffusion [Jj [21] . Despite numerous com- 
puter simulations the fundamental physics of hydrated electron diffusion still 
remains poorly understood [351 H31 121] • Grotthus type mechanisms, ad iabatic 
long-range hopping, and tunnelling effects have all been discounted. Instead, 
a polaron-like picture of the dynamics has emerged in which the high mobility 
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has been attributed to the instantaneous adiabatic response of the electron to 
fast solvent motions [22 | 123 ] [25 j [26] . 

In this work we demonstrate that the mechanism and kinetics of the diffusion 
process are well described by a transfer diffusion model [21 , 28J. We consider a 
transfer process involving the exchange of an extra-molecular electron between 
identical solvent molecules A x and A 2 , 

[M,^] + A 2 ^A 1 + [A 2 ,e m ]. (1) 

For a system populated entirely of A\ and A 2 migration of the excess charge 
occurs either by the above charge-transfer process and/or by the conventional 
diffusion of [A\, e~ q ] or [A 2 , e~ q ]. If the timescales associated with the diffusion of 
A and the electron-transfer process are similar in magnitude both processes may 
contribute significantly to the diffusion of charge. For example, voltammetric 
studies of electron conduction in inorganic solutions indicate that the charge- 
transfer process provides a significant contribution to the overall migration of 
charge [25] . 

We employ an adiabatic quantum-classical molecular dynamics method to 
investigate hydrated electron diffusion over the temperature range 298 K to 
400 K. In common with other works [22. 23, 24] the system is evolved within the 
Born-Oppenheimer approximation with the solvent molecules treated classically. 
The electron, confined to the ground state, is represented by a wavefunction that 
corresponds to the instantaneous nuclear configuration. Our methodology has 
been previously employed for the hydrated electron at different thermodynamic 
states [3ni[3T], in confined media [35], and in electron-cation pairs [3"3l l3"4l 1531 13"6"] . 
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2 Method 

Mixed quantum-classical molecular dynamics is used to simulate the excess elec- 
tron in bulk water [30l GZ] . The problem of the hydrated electron is reduced to a 
localised quantum mechanical problem within a classically evolving system [38) . 
Water- water interactions are treated classically, using the SPC model [39], with 
the solvent coordinate, S, evolved according to Newton's equation of motion. 
The electron-solvent interaction potential of a given configuration is termed 
Vint(r, S), where r denotes the electronic coordinates. The energies (E n (S)) 
and wavefunctions (^> n (r, S)) of the excess electron may be computed by solv- 
ing the time-independent Schrodinger equation, 

f e + \\ a t (r , S)] ip n (r , S) = E n (S) V>„ (r , S) , (2) 

where T c is the electronic kinetic energy operator. The electron-water inter- 
action is described by the pseudopotential of Turi and Borgis [ID] which has 
been designed specifically for liquid simulations. The principal advantages of 
this pseudopotential are its simple analytical form, and its good reproduction of 
both the equilibrium ground state energy and experimentally measured quanti- 
ties (e.g. optical absorption spectrum). 

The Born-Oppcnhcimcr approximation, from which we assume the separa- 
tion of timescales between the electron and the classical subsystem, allows for 
an adiabatic simulation of the excess electron. Thus the electron wavefunction 
is confined to the ground electronic state corresponding to the instantaneous 
configuration of the nuclear (solvent) coordinates. The quantum mechanical 
contribution to the forces acting upon the solvent from the excess electron are 
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computed according to the Hcllmann-Feynman theorem [5TJ 22] : 




d_ 
OS 




dr |^ n (r,S)| 2 — V- nt (r,S). 



(3) 



We compute only the contributions arising from the ground state (n = 0) of 
the excess electron. We note that the energetic difference between the ground 
and the first excited state of the excess electron is sufficiently large so as to ren- 
der any thermal excitation negligible. The adiabatic mixed quantum-classical 
methodology outlined above is sufficiently inexpensive so as to allow the rela- 
tively long simulation times (~0.5 nanoseconds) necessary for investigating the 
diffusive dynamics. Indeed one electron models, of the kind presented here, have 
been shown to provide a good approximation to multi-electron approaches, de- 
spite omitting the spreading of the excess electron density into the anti-bonding 
orbitals of hydrating water molecules [43;. For a more detailed presentation of 
our methodology the reader is referred to the relevant literature [30l |37j [40] . 

The excess electron was simulated in a bulk system (density 1.00 g/cm 3 ) of 
300 water molecules at 298 K. Simulations were performed within the canonical 
ensemble using the Nose-Hoover thermostat [44] . Periodic boundary conditions 
were employed and long-range forces treated using the Ewald summation tech- 
nique |45j . The Gear predictor-corrector algorithm was used to integrate the 
translational and rotational equations of motion [45, 46 . After equilibration the 
system was simulated for 500 picoseconds, with a timestep of 0.5 femtoseconds. 
Configurations were sampled every 5 femtoseconds. This process was repeated 
at progressively higher temperatures: 325, 350, 375, and 400 K. All simulations 
were performed at the same density of 1.00 g/cm 3 . Additional simulations were 
performed on bulk water systems without the excess electron. 



5 



3 Results and Discussion 

As previously demonstrated [SOJ SO], the hydrated electron is found to be a 
well-localized entity in agreement with experiment [471 148) and other simulation 
works [2H [231 [25J [571 EH] • As such the electron position may be taken as its cen- 
tre of mass, r c — (ip \r \ ip). Though the electron wavefunction is not propagated 
in time we may examine the electron's dynamical behaviour semi-classically 
through its centre of mass. Figure [T] compares the mean-squared displacement, 
(Ar c (£) 2 ), of the hydrated electron with bulk SPC water at 298 K. The propor- 
tionality between (Ar e (i) 2 ) and t after 1 ps, indicates that the electron under- 
goes Brownian-like diffusion. Assuming Brownian motion, the related diffusion 
coefficient is given by the Einstein relation, 

d(Ar c (t) 2 ) 

-A_pi = 6Dt , (4 ) 

where Ar e (t) is the total displacement between t = and t, and D is the respec- 
tive coefficient of diffusion. The electron diffusion, found here to be 0.60 A 2 ps _1 , 
is approximately 1.6 times faster than the SPC solvent, 0.38 A 2 ps _1 . In con- 
trast to early simulation work the electron diffusion is clearly greater than the 
solvent [531 [23]. Using a similar methodology with the rigid, polarisable TIP4P 
model for water [50], Staib and Borgis have estimated the electron diffusion at 
0.4 A 2 ps _1 [57] . Both here and in the work of Staib and Borgis fast diffusion 
occurs without the inclusion of tunneling effects. Long-range hopping effects are 
also absent. Conductometric experiments estimate the hydrated electron diffu- 
sion at 0.49 A 2 [TJdT] approximately 2.1 times faster than water self-diffusion, 
0.23 A^s" 1 [51]. 

A direct examination of the electron position reveals not only the long-time 
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diffusive dynamics but aiso the shorter-time, osciliatory motions (Figure |2p,) . 
The instantaneous response of the electron to the solvent coordinate correlates 
the electron and solvent dynamics. As such we may assume that the timescales 
of the electron dynamics mirror the individual and collective modes of the sur- 
rounding water molecules. We determine the frequencies of the electron motion 
from the spectral density, S v (u>), of the electron centre of mass velocity auto- 
correlation function C v (t), 



The electron velocity, v(£), is computed from the particle displacement between 
sample intervals. S v (ui) is shown in Figure for the hydrated electron at 
298 K. The peak frequencies 376 cm" 1 and 653 cm -1 (extracted through cubic 
interpolation) lie within the region characterised by intermolecular vibrations 
in bulk water. Low frequency Raman spectra of the hydrated electron exhibits 
two main peaks at around 450 cm -1 and 700 cm -1 [S"2"j . 

The radial distribution function between the electron position and surround- 
ing water molecules reveals a hydration shell comprising approximately 4 water 
molecules |40j . This hydration shell is defined according to the first minimum 
(at approximately 3.7 A [3UJH0]) in the e~-water RDF. To further investigate 
the origin of short-time electron dynamics we have evaluated the librational fre- 
quencies of both the hydrating and 'bulk' water molecules. Figure [3] shows the 
librational frequencies of hydrating (r c _o < 3.7 A), and 'bulk' (r _o > 3.7 A) 
water molecules at 298 K, computed according to the methodology of Bopp [S3] 
(see appendix). The peak frequencies, extracted from the interpolated data, 




(5) 



C v (t) 



(v(O)-vW) 
(v(O)-v(O)}' 



v(t) 



r e (t + At/2) - r e (t - At/2) 
At 



(6) 
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are given in Table [TJ The dynamical behaviour of water molecules characterised 
as 'bulk' is not greatly perturbed by the excess electron. For bulk water the 
restricted rotations (or rocking motio ns) that form the librational motions are 
manifested in vibrational spectra by the two bands, L\ and £2- Estimates of 
peak frequencies from IR and Raman spectra vary; with L\ located around 
380 cm -1 and L 2 around 670 cm -1 at ambient temperature [541 155] . R x and 
R z constitute the lower frequency L\ band and R y corresponds to the broader 
L2 band. 

The significant red-shift in the librational frequencies from 'bulk' (Figure [3]d) 
to hydrating water (Figure [3^,) indicates that the electron perturbs the libra- 
tional dynamics of the hydrating water molecules. This red-shift has been ob- 
served in both experiments and simulations [33]. A blue-shift in librational 
frequencies has been observed for hydrated ion systems [56l |57l [58] , where the 
magnitude and direction of the frequency shift indicate ion-water interactions 
greater than water-water interactions. In this instance the weaker electron- water 
interaction red-shifts the librational frequencies. 

On comparing the computed librational frequencies with the electron spec- 
tra in Figure [2] the lower frequency peak (at 376 cm^ 1 ) can be associated with 
hydrating water molecules (mainly R x and R z librations) whilst the higher fre- 
quency peak (at 653 cm -1 ) can be attributed to non- hydrating water molecules 
(mainly the R y libration). These results confirm the long held notion, from ex- 
perimental H31 031 [59] and simulation work [T5], that the electron dynamics is 
dominated by water molecule librational motions. Furthermore, the direct cor- 
respondence between the librational and electron frequencies shows the electron 
to be highly responsive to individual water molecule motions. As such electron 
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diffusion need not necessarily require collective water dynamics. 

Whilst the hydrated electron is a well-localised entity the fast oscillatory dy- 
namics of the electron position suggest a highly fluxional local potential energy 
surface. In accordance with the conventional structural picture of electron hy- 
dration [47] this potential energy surface is expected to be characterised by local 
minima near the oriented dangling H atoms of hydrating water molecules [40j . 
The short-time dynamics observed in the particle coordinate (Figure [2]) may 
then be understood as oscillatory movements between dangling H atoms. Elec- 
tron motion being prompted by the fluctuating electrostatic field of the solvent 
librations. These rapid displacements between H atoms would be characterised 
by a free energy barrier on the order of k^T. The barrier itself being a func- 
tion of the solvent coordinate (both the position and orientation of surrounding 
water molecules). 

Computer simulations of Park et al. have demonstrated the importance of 
dangling hydrogen atoms in hydrated electron dynamics |26| . They describe a 
process in which the electron, driven by dangling H atoms, responds to thermal 
fluctuations in the H-bond network. The destruction and formation of dan- 
gling H atoms in the surrounding H-bond network pushes and pulls the electron 
between water molecules. We assume that the diffusive process consists of a 
sequence of displacements between dangling H atoms. As such, the diffusion 
may be described by the exchange process identified in equation [I] To compute 
the decay of [A±, e~ q ] in equation [l] we define the following simple residence-time 
correlation function, 

CR{t) ~ (P(o)P(o)>- (7) 

The operator P(t) is 1 if the electron, initially attached to a given H atom at 



t = 0, is still attached at a later time t. The electron is defined as 'attached' 
to the closest H atom. Cn(t) is thus the residence time of the electron on a 
given H atom. The decay of Cji(t) is thus the rate at which the electron moves 
between dangling H atoms. Cn(t) is shown in Figure|4]for 298 K. The initial fast 
decay of Ca(t) (t < 0.05 ps) reflects the fast oscillatory short-time dynamics. 
At longer times (t > 0.2 ps) the short-time dynamics averages out revealing 
the slower long-time diffusive dynamics. Also shown in Figure [4] is r c _H, the 
average distance between the electron and the H atom to which it was initially 
attached. From Figure [4] we may estimate, by visual inspection, that electron 
displacements between dangling H atoms are within the range 0.7-1.0 A. 

The long-time proportionality of 1/Cn(t) with time (Figure [5]) confirms the 
second-order rate law typical of self-exchange processes. The gradient at long 
times gives the phenomcnological rate constant, fc cxc , associated with the ex- 
change process. The values of fc exc are listed in Table [III Noticeably, fc oxc appears 
to be faster than the 1 ps timescale often associated with hydrogen-bond rear- 
rangements in bulk water (6(2 EH [62] . 

The calculated electron diffusion, D c \ cc , and the exchange rate constant, 
k cxc , are plotted in Figure [6^,. Figure [6]d shows that our model reproduces the 
Arrhenius-like electron diffusion. The Arrhenius behaviour of -D c iecj an d k exc 
constrasts strongly with the non- Arrhenius behaviour of the water self-diffusion. 
Clearly, neither the exchange process nor the electron diffusion are driven by 
the water diffusion. The (constant volume) activation energies for the electron 
diffusion and the charge transfer process are 8.9±0.7 kJ mol" 1 and 10.0±0.4 kJ 
mol -1 respectively. The close agreement strongly suggests that the electron 
diffusion is well described by the exchange process in equation [T] Experimental 
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estimates of the activation energy are much larger at around 20.0 kJ mol -1 [7]. 

As a further test we assume proportionality between fc oxc and D c \ cc . We 
consider the Brownian behaviour of the electron as a simple random walk in 
which the step frequency is given by fc eX c! 

e 2 fc CX c 

-Dcaic = — -p, — • (8) 
o 

The parameter e is the electron displacement accompanying the exchange pro- 
cess and is extracted from a linear fit of k exc against I? i cc (Figure 0. The 
displacement, e, is found to be 0.85 A and is in good agreement with our earlier 
estimates from r c _n (see Figure [4]). The good agreement between D ca i c and -D e iec 
(see Table 2) confirms that k exc is indeed the phenomenological rate constant 
describing the hydrated electron diffusion. 

In agreement with earlier works on electron diffusion [2TJ [2U [25] these results 
confirm the central importance of both the fast reorientational solvent motions 
and the instantaneous electron response. They also demonstrate that electron 
diffusion occurs rapidly at ambient temperatures purely under the influence of 
thermal fluctuations, and as such tunnelling processes arc unnecessary to explain 
its high mobility [23] . The model presented here reproduces the Arrhenius-like 
behaviour of the diffusion but underestimates the activation energy. The experi- 
mental Raman spectra of the hydrated electron is also qualitatively reproduced. 
A quantitative description of the electron dynamics requires a more accurate de- 
scription of the fast solvent dynamics, including both intra-molecular degrees of 
freedom and polarisability. Comparison of computed and experimental spectra 
provides a clear test of the strength of the model [431 E2] ■ 

The clear importance of librational dynamics on the electron diffusion raises 
an apparent paradox. The fast diffusion of the electron is enabled by the libra- 
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tional dynamics of the solvent. As such it may be expected that the increase in 
diffusion with temperature is due to an increase in the frequency of librational 
motions. However, librational frequencies decrease with temperature |55) . The 
increase in diffusion appears to result from greater librational amplitudes and a 
more disordered hydration layer. Both may contribute to the fluxional character 
of the free energy barrier associated with electron displacements. Crucially, the 
respective back reaction of equation [l] is expected to diminish at higher tem- 
peratures as orientational configurations within the coordination sphere change 
much more quickly. 

4 Conclusion 

We have employed mixed quantum-classical molecular dynamics simulations to 
investigate the short and long-time dynamics of the hydrated electron within 
the temperature range 298-400 K. The short-time dynamics is characterised by 
oscillatory motions between dangling H atoms. This is a thermally activated 
process promoted by the librational dynamics of surrounding water molecules. 
Displacements between dangling H atoms constitute the fundamental step in the 
electron diffusion. We have demonstrated that the associated mechanism and 
kinetics can be understood as a transfer diffusion process in which the electron 
is exchanged between neighboring water molecules. The Brownian diffusion 
may then be modelled as a simple random walk in which the step frequency 
is given by the rate constant for the exchange process. The results presented 
here provide a consistent picture of both the short and long-time dynamics of 
the hydrated electron. They facilitate further investigations into the a pplica- 
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bility of traditional ion diffusion models, water reorientation in the presence 
of anions [63]. and the seemingly anomalous hydrated electron diffusion at low- 
temperatures [7] and within confined systems |32j . 
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6 Appendix 

The three librational motions are labelled R x , R y , and R z . The librations are 
approximated as vibrational modes corresponding to the instantaneous rotations 
of the water molecules about their three principal axes. Using the methodology 
of Bopp [S3], we define two sets of unit vectors to describe the vibrational 
coordinates from which the modes are constructed. The first set of unit vectors, 
Vi and v 2 , assigned to H atoms 1 and 2 respectively, are perpendicular to the 
O-H bonds, within the plane of the molecule, pointing outwards. The second set 
of vectors, pi and P2, are defined perpendicular to the plane of the molecule. 
The projection of the H atom velocities along the unit vectors v and p are 
labelled V\, V2, Pi, and P 2 - We construct the normal coordinates R x , R y and 
R z for a water molecule lying within the yz-plane wit h its C2 axis aligned to 
the z-axis, 

R* = Vi- V 2 

Ry =Pl+P 2 (9) 

Rz=Pl-P2 ■ 
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R x decribes rotation about the x axis; it is the in-plane swinging motion of 
the water molecule, often termed 'rocking'. R y , commonly referred to as the 
'wagging' motion, is a rotation about the y-axis during which the H atoms swing 
backwards and forwards across the molecular plane. Finally, R z describes the 
twisting of the water molecule about the symmetry axis. The spectral densities 
derived from the autocorrelation of the above normal coordinates provide the 
respective librational spectra. For further details readers arc referred to the 
relevant literature j53j [56] . 
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Hydrating (cm 1 ) 


Bulk (cm" 1 ) 


Bulk Experiment (cm 1 ) 


Rx 


379 


462 


380 


Ry 


411 


656 


670 


R z 


325 


450 


380 



Table I: Computed librational frequencies for hydrating (r c _o < 3.7 A) and 
'bulk' (r c _o > 3.7 A) water molecules. Experimental values for bulk water are 
approximate |55j . Peak frequencies are extracted from a cubic spline interpola- 
tion of the data set. 
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Temperature 


Dh 2 o 




^exc 


e 




(K) 




(A^s- 1 ) 


(ps- 1 ) 


(A) 


(A^s- 1 ) 


298 


0.39 


0.60 


5.00 


0.85 


0.60 


325 


0.57 


0.94 


6.55 


0.85 


0.79 


350 


0.74 


1.12 


9.19 


0.85 


1.11 


375 


0.87 


1.35 


11.09 


0.85 


1.34 


400 


0.97 


1.58 


13.85 


0.85 


1.67 



Table II: The self-diffusion coefficient of SPC water, -Dh 2 o, the electron diffu- 
sion coefficient, D e \ ec , the electron exchange rate constant, fc cxc , the estimated 
random walk step distance, e, and the random walk model diffusion, D ca i c . As 
Cr(*) is dimensionless, the second order rate constant, fc exc , is given as ps -1 . 
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Figure 1: Mean-square displacement of hydrated electron centre-of-mass com- 
pared with that of bulk SPC water. 
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Figure 2: Top, x-coordinate of the electron centre of mass as a function of 
time. Bottom, power spectrum of the electron velocity autocorrelation calcu- 
lated according to equations [5] and [6j The curve is obtained by cublic spline 
interpolation. 
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Figure 3: Librational spectra of hydrating, (top, r _o < 3.7 A), and 'bulk' 
(bottom, r e _o > 3.7 A) water molecules. See appendix for computational details. 
Curves are obtained by a cubic spline interpolation. 
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Figure 4: Residence time correlation function, Cu{t), computed at 298 K. Also 
shown is r c _n, the average distance between electron and H atom during the 
electron displacement. 
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Figure 5: Plot of l/C R (t) against time at 298 K, 350 K and 400 K. Straight 
lines are linear fits to the long-time portion of l/Cn(t). 
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Figure 6: Top, a plot of the electron diffusion coefficient, D e \ ec , and the exchange 
rate constant, fc oxc , with temperature. Bottom, Arrhenius plots for D elcc , fc cxc , 
and the bulk diffusion of SPC water, -Dh 2 o- 
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k e XC (P S "') 

Figure 7: A plot of -D i cc against fc oxc . Also shown is the line of best fit through 
the origin. 
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